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INTRODUCTION 

This  report  studies  microwave  backscattering  from  a  for¬ 
est  canopy.  The  canopy  is  modeled  by  lossy  dielectric  discs 
having  random  orientation.  Explicit  expressions  are  derived 
for  the  horizontal,  vertical  and  cross  polarized  backscatter¬ 
ing  coefficients.  The  results  are  then  compared  with  experi¬ 
ment. 

This  work  was  motivated  by  the  need  to  relate  radar  re¬ 
turn  to  the  characteristics  of  the  scattering  objects.  The 
situation  is  particularly  complex  in  the  case  of  vegetated 
terrain  which  consists  of  an  ensemble  of  many  highly  irregu¬ 
lar  objects  placed  in  a  more  or  less  random  fashion.  Using 
the  scattered  return  from  a  vegetated  terrain  one  would  like 
to  obtain  such  information  as  height,  density,  average  leaf 
and  branch  orientation  so  that  the  vegetation  can  be  classi¬ 
fied.  This  type  of  information  is  important  for  military 
analysis  of  the  terrain.  In  addition,  if  he  terrain  has 
been  determined  adequately  from  an  electromagnetic  point  of 
view,  the  backscatter  information  can  be  used  to  predict  ra¬ 
dar  returns  at  other  angles  and  frequencies  not  measured. 

This  would  be  of  great  aid  to  the  development  of  manageable 
data  bases  for  radar  simulations. 

The  above  applications  have  served  as  a  motivation  for 
the  development  of  electromagnetic  models  for  the  vegetated 
terrain.  These  models  have  been  constructed  by  replacing  the 
vegetated  region  with  a  random  medium  whose  statistical  char¬ 
acteristics  are  related  to  the  physical  quantities  of  the 
medium.  The  models  can  be  divided  into  two  categories: 


continuous  and  discrete.  In  the  continuous  case,  the  vegeta¬ 
tion  is  modeled  by  assuming  that  its  permittivity  e (x)  is  a 
random  process  whose  moments  are  known.  The  average  back- 
scattering  cross  section  is  then  calculated  from  a  knowledge 
of  the  statistics  of  e (x) .  Usually  it  is  only  the  mean  and 
correlation  of  the  permittivity  that  are  required.  The  anal¬ 
ysis  of  the  problem  proceeds  by  first  calculating  the  mean 
wave  in  the  random  medium  [Keller,  1962]  and  [Tatarskii  and 
Gertsenshtein,  1963] .  The  mean  wave  allows  the  characteriza¬ 
tion  of  the  vegetation  by  an  equivalent  dielectric  constant. 

The  scattered  field  is  then  calculated  by  using  single  scat¬ 
tering  theory  in  this  equivalent  medium  -  a  technique  employ¬ 
ed  by  Rosenbaum  and  Bowles  [1974]  .  Particular  applicaiton  of 
the  method  to  vegetated  medium  has  been  made  by  Lang  [1974],  Hevenor 
[1976] ,  Fung  and  Fung  [1977] ,  Fung  and  Ulaby  [1978] ,  Fung  [1979] ,  and 
Zuniga,  et  al,.  [1979].  The  theory  is  scalar  in  origin  and 
thus  does  not  provide  cross  polarization  information.  In  ad¬ 
dition,  a  correlation  function  for  medium  fluctuations  must 
be  assumed.  As  yet  the  detailed  structure  of  this  correla¬ 
tion  function  has  not  been  related  to  medium  properties. 

In  this  report,  we  have  adopted  the  alternative  approach 
-  modeling  by  discrete  random  media  techniques.  Here,  the 
individual  objects  -  such  as  leaves  -  are  characterized  by 
their  scattering  cross  sections  or  dipole  moments.  Each  ob¬ 
ject  is  then  given  a  random  placement  and  orientation.  The 
analysis  proceeds  in  a  similar  manner  to  the  analysis  of  the 

t 

continuous  case.  First  the  mean  wave  in  the  medium  is  found 
by  using  techniques  developed  by  Foldy  [1945],  Lax  [1951], 


Twersky  [1962] ,  and  Keller  [1964]  .  Then  an  equivalent  die¬ 
lectric  tensor  is  found  and  the  backscattering  is  computed  by 
using  single  scattering  in  the  equivalent  medium.  In  appli¬ 
cations  for  vegetation  modeling,  the  method  has  been  used 
previously  by  Lang  and  Sokolakis  [1979]  for  spherical  scattered s.  We  new 
extend  the  method  to  arbitrary  shaped  scatterers.  Now  we 
find  that  the  equivalent  medium  is  anisotropic  in  nature. 

This  leads  to  zeroth  order  depolarization  of  the  incident 
wave.  In  addition,  the  backscattering  coefficients  are  ex¬ 
plicitly  related  to  scatterer  volume  and  dielectric  constant. 

Previous  work  using  discrete  scatterers  to  model  vegeta¬ 
tion  has  been  done  by  Du  and  Peake  [1969] ,  but  they  employed 
single  scattering  (Born  approximation)  without  introducing  an 
equivalent  medium.  Thus  they  did  not  take  into  account  the 
decay  of  the  incident  wave  in  the  vegetation.  This  limits 
their  theory  to  much  lower  frequencies  and  thin  layers  of 
vegetation . 


PROBLEM  FORMULATION 


Consider  the  problem  of  scattering  of  time  harmonic 
electromagnetic  waves  from  N  discrete  scatterers  located  in 
a  volume  V  as  is  shown  in  Figure  1.  The  particles  are  all¬ 
identical  and  each  has  volume  V  ,  relative  dielectric  constant 

? 

er  and  free  space  permeability  .  It  is  assumed  that  the 
background  medium  is  free  space  having  permittivity  and 
permeability 

.  tii 

The  position  of  the  i  particle  is  specified  by  the 
vector  extending  from  an  origin  0  to  the  center  of  that 
particle.  The  particle’s  center  is  located  by  the  center  of 
the  smallest  circumscribed  sphere  in  which  the  particle  can 
be  placed.  Although  the  particles  are  identical  they  have  a 
rotation  with  respect  to  a  fixed  direction.  The  rotation  for 
the  i particle  is  specified  by  =  ,  <*k)  where  0^  and 

are  polar  and  azimuth  angles  respectively  with  O<_0^<tt  and 
0  <  $  .  <  2  7T  . 

—  i— 

The  electric  field  obeys  the  vector  wave  equation 

VxVxE  -  k^er(x)E  -  ioy0  J  (1) 

where  a  time  dependence  e  IUJt  has  been  assumed.  In  (1)  = 

u)  / r.  q  Uq  is  the  free  space  wavenumber  and  J  is  the  current 
density  of  the  source.  The  relative  dielectric  constant 
e^(x)  can  be  expressed  in  terms  of  individual  particles  by 
employing  translations  and  rotations  of  the  particle  located 
at  the  origin.  Let  us  assume  that  a  particle  located  at  the 
origin  is  characterized  by  the  function  U(x)  where 


anriwfpwr j^fwnuT) 


U  ( x)  = 


1  ,  xeV 

~  I 

0  ,  x^Vy 


(2) 


Using  (2)  we  express  er(x)  as 


N 

e,.(x)  =  1  +  A  l  U(x-X. ,ft.)  , 

r  -  ~  -i  ”i 


A  =  er-l  (3) 


where 

U(x#ft)  =  U  (B(ft)  -x)  (4) 

Here  U(x,ft)  is  the  function  U (x)  rotated  by  ft  and  R'ft)  is  a 
rotation  dyadic. 

We  will  find  it  convenient  to  express  (1)  and  (3)  in  a 
more  abstract  notation.  We  have 


where 


N 

(L  -  Z  V.  )  *E  -  a 
i=l  ~ 


L  *  VxV:cI  -  kQI 


V. 

=1. 


kQAU(x-Xj  ,ft^)  I  , 


3.  = 


iw^0- 


(5) 

(6) 

(7) 


Here  I  is  the  unit  dyadic  and  a  can  be  viewed  as  a  normaliz¬ 
ed  source  term.  At  times  it  will  be  convenient  to  write _ 

E(x)  *  Eq(x)  +  Eg(x)  (8) 

where  E^(x)  is  the  solution  to  (5)  when  no  scatterers  are 
present,  i . e.  , 

L*Eq (x)  =  a  (9) 


and  E  (x)  is  the  scattered  field  from  the  particles. 

— 3  — 


"9J 
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SINGLE  SCATTERED  -  TRANSITION  OPERATOR 

Before  considering  the  N  particle  scattering  problem,  we 
will  consider  scattering  from  one  particle  located  at  the 
origin.  Putting  N=1  in  (5)  with  and  we  have 

(L  -  V)  'e  =  g_  ,  V  ~  AkgO(xfQ)  (10) 

where 

e  =  eQ  t  es  ,  L-e_0  =  5  (11) 


and  e  is  outgoing  as  [x|-*°°.  We  have  used  the  small  e  nota- 
tion  for  the  field  here  to  remind  us  that  there  is  only  one 
scatterer  present. 

If  we  use  (11)  in  (10) ,  we  obtain 

=  ¥'£  ( !2 ) 


From  (12)  we  see  that  the  term  on  the  left,  V*e,  can  be  view¬ 
ed  as  the  source  of  the  scattered  field.  We  write 


9~  =  v’e  (13) 

where  is  an  equivalent  source  term.  Since  Y-0  when 

x^Vp,  the  sources,  5eg/  exist  inside  the  particle  boundaries. 

It  is  more  natural  to  think  of  the  equivalent  sources  as 
being  caused  by  the  incident  field  e..  Because  Maxwell's 


equations  are  linear,  we  can  write 

g  =  l’e. 

-^eq  -  —0 


■  (14) 


where  the  dyadic  operator  T  is  known  as  the  transition  opera¬ 
tor  in  the  scattering  literature  [Lax,  1951].  Now  using  (13) 
and  (14)  in  (12)  and  multiplying  through  by  L  1 ,  we  have 


eg  -"•  i”'L‘I*i0  (15) 

Thus  the  knowledge  of  I  completely  characterizes  the  scatter¬ 
ing  properties  of  the  particle.  The  operator  I  is  related  to 
the  dyadic  scattering  amplitude  of  the  particle  and  for  dipole 
scatterers  T  can  be  determined  from  the  polarizability  of  the 
particle.  Thus  T  is  directly  connected  with  quantities  of 
physical  interest. 

The  transition  operator  is  a  linear  bounded  operator  and, 
as  a  result,  can  be  expressed  in  integral  form: 

(x)  ~  T  •  e^  =  Jdx  '  t(x,x' )  •  e^(x')  (16) 

where  the  limits  for  the  integral  extend  over  all  space.  One 
can  show  that  £  is  0  when  x  and  x'  are  outside  the  particle 
[Frisch,  .1968],  i.e., 

t(x,x')  =  0  ,  x^Vp  or  xWp  (17) 

The  property  follows  directly  from  the  fact  that  the  equiva¬ 
lent  sources  for  the  scattered  field  are  located  within  the 
particle  boundaries. 

We  will  now  represent  %  in  terms  of  plane  waves.  The 
representation  of  £  can  be  directly  related  to  the  dyadic 
scattering  amplitude.  We  proceed  by  representing  e_c(x)  by 
its  Fourier  transform,  putting  this  in  (16)  and  taking  the 
Fourier  transform  of  (16) .  We  obtain 

2eg00  =/  dk't(k,k’)  *i0(k>)  _  (18) 

where 

t  (k  ,k  ' )  =  — --t?  /dxdx  '  t  (x,x' )  e*"1  (  — }  (19) 

(2tt)  J 

In  (18)  we  have  used  the  notation  that  h  is  the  Fourier 


i 


Stf”' - 


transform  of  h  •  More  specifically 

/+ 

h  (k)  =  J  dx  h(x)e  i-’~  (2  0) 

Inverting  (19)  the  transition  kernel  t  can  be  expressed  in 
terms  of  its  plane  wave  representation  t: 


t(x,x*)  =  — fdkdk'  t  (k,k:)el(~*-~-' (21) 

”  .  (  2  TT )  J  J 


The  dyadic  scattering  amplitude  will  now  oe  defined. 
Consider  a  plane  wave  incident  upon  a  scatterer  located  at 
the  origin.  An  arbitrary  plane  wave  can  be  decomposed  into 
two  mutually  orthogonal  linearly  polarized  plane  waves.  The 
polarization  directions  are  taken  as  a°  and  S_°  where  cx°  and 
j3°  are  orthogonal  unit  vectors  with  o°  and  6_°  being  perpen¬ 
dicular  to  the  direction  of  propagation.  The  two  incident 
waves  are 


“^0— *— 

en(x,i;q)  =  £°e  ,  qela,6}  (22) 


-Q  v 

where  i.  is  a  unit  vector  in  the  direction  of  incidence.  It 
is  more  convenient  to  consider  both  polarizations  simultan¬ 
eously  so  we  introduce  the  dyadic  incident  wave  [Twersky,  1967] 
—  ^  “  S-0  ^  a ^  9Lq  +  ( £/ L?  ^  §Lo  (23) 


“  (^0^)  +  -0^0}  6 


ikQ_i  'x 


(24) 


ik0i-2 


-  (I  -  ii_)  e  ”  (25) 

The  dyadic  scattered  field  from  the  particle  is  given  by 

es(x,j.)  =  £s(x,i:a)a0  +  £s (x ; S ) (26) 


where  e  (x,i;q)  is  the  scattered  field  due  to  polarization  q. 


The  dyadic  scattering  amplitude,  f,  is  defined  in  terms  of  che 
asymptotic  expression  for  e_  in  the  radiation  zone.  We  have 


r 


X  +00 


(21) 


ik0  !  x  I 

=s  *—  ±}  *  ^7x] 

where  0_  is  a  unit  vector  in  the  x  direction,  £=x/|x|* 

The  relationship  between  f  and  t  can  be  found  by  employ¬ 
ing  (15)  for  large  |xj  (Appendix  A) .  The  result  is 

|(0,i)  -  2tt2(J-0  0)  -|‘(k Q0,k  i)  •  (J-i  i)  (28) 


From  this  relation,  we  see 

0*|  =  0,  f * i  *  0  (29) 

Thus  f  is  a  four  component  tensor  ~  all  combinations  of  two 
incident  polarizations  and  two  scattered  polarizations.  We 
also  note  that  £  does  not  completely  determine  t  but  only 
partially  specifies  it.  In  particular,  a  knowledge  of  f  for 
a  free  space  wave  number  kQ  only  determines  £(k,k')  at  |k|  = 

|k' |=kQ;  also  only  four  of  the  nine  components  of  t  in  the 
polarization  directions  are  determined. 

Before  concluding  this  section,  transition  operators  for 
particles  not  located  at  the  origin  will  be  needed.  As  before, 
the  equivalent  sources  for  a  particle  located  at  can 

be  related  to  the  incident  field .  It  follows  that 

^(.x)  =T  i*eQ  =  ft.  _(x,x')  -^(x'Jdx'  (30) 

By  shifting  the  sources  and  the  incident  field  to  the  origin 
t.  can  be  related  to  t.  One  finds 

=i  — 

t. (x,x’)  *  £{x-X. ,x'-x. )  .  (31) 

KX  .  —  — X  X 

Note  that  throughout  the  discussion  the  dependence  of  t  on 
rotations  has  been  suppressed  for  convenience  of  notation. 
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COHERENT  FIELD 

In  this  section  we  will  develop  an  approximate  equation 
for  the  coherent  field  by  employing  the  Foldy  approximation 
[Foldy,  1945] .  The  equation  is  in  terms  of  the  transition 
operator  and  thus,  when  the  scattering  amplitude  is  known, 
the  equation  is  completely  specified.  After  the  equation  has 
been  derived  it  is  pointed  out  that  outside  V  the  coherent 
fields  obeys  Maxwells '  equations  with  free  space  permittivity 
and  permeability.  Inside  V,  the  coherent  fields  obey  Maxwells 
equations  with  free  space  permeability  and  a  macroscopic  per¬ 
mittivity  that  is  inhomogeneous,  anisotropic  and  spatially 
dispersive . 

Before  discussing  the  coherent  field,  the  statistics 
that  govern  the  particles  position  and  rotation  must  be  specified.  It 
will  be  assumed  that  the  position  vectors  X^,  i=l...N  and  ro¬ 
tation  vectors  £2^,  i=l.  .N  are  random  variables  that  are  spe¬ 
cified  by  a  5n  dimensional  distribution  function.  In  addi¬ 
tion,  it  .is  assumed  that  interchanging  particles  leaves  the 
distribution  function  unaffected.  From  this  general  distri¬ 
bution  function  we  can  obtain  the  probability  density  func- 
th 

tion  for  the  i  particle.  It  is 

Px  s  i=l - N  (32) 

— i — i  — 

where  tu  =  ( 8 ,  4> )  .  In  (32)  we  have  explicitly  noted  the  fact 
that  the  particles  are  identically  distributed  by  omitting 
the  index  i.  on  the  left  hand  side  of  (32)  .  We  will  assume 
that  the  particles  location  and  rotation  are  independent,  thus 

Pxq(x^)  =  Px(x)pn(o;)  (33) 


mmmvm+mi  imwpuwwpmwm  i 
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with  the  usual  property 


ypv (x) ax 
V  $  ~  ~ 


=  1 


f. 


The  particle  density  is  defined  by 


p(x)  =  N  p  (x) 


p0(w'dii  =1  (34) 


(35) 


so  that 


/. 


p (x) dx  =  N  . 


(36) 


In  addition  to  the  one  particle  density  -  when  treating  the 
correlation  of  the  field  -  the  two  particle  density  will  be 
required.  We  have 

A  A  A  A  A  A 

Px.x.ft.$2.  PXXQfl(-,-,-'-)  =  (37) 

-  —  — 

i , j  «  1 ,  . .  .N 

In  (37)  we  have  assumed  that  the  i  and  j  particles  are 
independent.  The  independence  assumption  is  valid  when  the 
particles  are  sparsely  distributed;  the  case  we  intend  to 
treat . 

We  will  now  develop  the  approximate  equation  for  the 
coherent  field.  We  start  by  noting  that  the  total  field  E 
can  be  thought  of  as  a  sum  of  the  incident  field  EQ  plus  a 
sum  of  the  fields  scattered  from  each  particle,  .  We 

have 

N  .  ,  .  . 

E  =  E  +  L  E*1'  (38) 

i=l  S 

The  total  field  incident  on  the  i^  particle  is  called  the 
effective  field  and  is  denoted  by  .  Thus  repre¬ 

sents  the  equivalent  sources  generated  by  the  incident  field 
in  the  it^‘  particle  and 


E(i) 
—  s 


J"1.Ti-E(i) 


39) 


Using  (39)  in  (38),  we  have 


N 


-1  ( i) 

E  =  EQ  -4-  Z  L  ~  •  T  .  •  E 

-  ■->  i = 1 


(40) 


This  is  the  equation  that  we  wished  to  obtain. 

Now  we  average  this  equation.  The  result  is 


N 


<E>  =  EA  +  Z  L_1  ■<  T  .  *E(l) 
"  ~°  i-1  =  =1  - 


(41) 


To  obtain  an  approximate  equation  for  the  mean  we  follow 
(Folay,  1945]  and  assume 


E  ^  -  <E> 


(42) 


(i) 


This  means  that  the  random  quantity  E  is  to  first  order 
equal  to  a  deterministic  quantity,  i.e.,  to  first  order  it 


is  a  ergodic  quantity.  Using  (42)  in  (41)  and  noting  that 
(i) 


<TV*E  >-<T^  *<  E>>=<  T^>  *<  E>  we  have  the  approximate;  equation 
for  the  mean  field 


N 


-1 


<E>  =  En  +  Z  L  *<T .  >  <E> 
u  i=1  =  -x 


(4  3) 


Denoting  explicitly  the  dependence  of  TU  upon  and  ,  av¬ 
eraging  and  then  using  (3  3)  ,  we  have 


T.>  =  <T(x.#n.)>  =  /  ds  /  dw  pvnls,u)I(s, 
-1  -  -1  -1  yv  -y4lT  -  XU 


w) 


=  /  ds  pv(s)T(.s) 

*V  ~  ~  ~ 


(44) 


where 


I(S)  =  f  du 


4  Tt 


Pq  (i£)  T  (s  ,w) 


(45) 


In  (45)  the  bar  over  T  has  been  used  to  indicate  an  average  of 


la  I  -J|  jM 


angular  varies  only.  By  putting  (44)  in  (43) ,  by  noting  that 
the  cattered  terms  are  identical  and  by  introducing  p (s)  via 
( 35 ) ,  we  obtain 

<£>  =  E.  +  /  ds  p  (s)L_1-T(s)  -<E>  (46) 

Multiplying  from  the  left  by  ^  and  using  (9) ,  we  get 

«t-<E>=  g  (47) 


where 


ds  p(s_)T(£) 


(48) 


This  is  the  equation  for  the  coherent  field. 

The  arguments  that  have  led  to  the  approximate  equation 
(47)  have  been  largely  heuristic.  The  essential  approximation 
is  contained  in  (42)  where  the  effective  field  is  assumed' 
approximately  equal  to  the  average  field.  Although  we  will 
not  discuss  the  conditions  under  which  (42)  is  valid,  it  will 
be  shown  elsewhere  that  the  approximation  is  valid  when  the 
fraction  of  volume  occupied  by  the  particles  is  small  com¬ 
pared  to  the  total  volume,  i.e.,  NV ^/V<<1.  We  shall  refer  to 
a  distribution  of  scatterers  satisfying  this  condition  as  a 
sparse  distribution. 

Before  proceeding  we  will  write  the  equation  for  the  mean 
in  more  concrete  form.  Using  (30)  and  (31)  in  (47)  and  (48) , 
we  obtain 

£,-<E(x)  >  -  J'  ds*  /*dx‘  P  (s)  t  (x-s,  x  '  ~s)  *<  E  (x 1  )  >  -  g  (49) 

V  J 


where 


Here  the  kernel  t(x,x';u)  is  the  same  as  given  in  (31) ,  how¬ 
ever  we  have  explicitly  shown  its  dependence  on  the  angular 
coordinate  t±>. 

We  can  now  use  (49)  to  obtain  a  macroscopic  form  of 
Maxwell's  equations.  First  averaging  the  Faradary 's  law 
equation;  we  have 

Vx<E(x)>  =  itqjQ<H(x)>  .  (51) 


Then  by  using  (51)  in  (49),  we  obtain  the  macroscopic  Ampere's 
law  equation. 

Vx<H(x)>  =  J  -  ico<D>  ,  <D>  =  £q|-<e>  (52) 

when  e  is  a  macroscopic  permittivity  operator  which  describes 
the  average  behavior  of  the  medium  and  J=g/ ( iojpQ )  •  It  is 

e  =  I  +  i-*-  /  ds  /dx*  p  ( s)  t  (x-s  ,  x  '  ~s )  .  (5  3) 

=  ~  *o  v  J  ~ 

This  expression  simplifies  to  |  (free  space)  when  xjz'V.  To 
see  this  we  note  that  when  x/V ,  we  have  x-s£V  since  sr V . 

Now  using  (17)  we  have  t=0 .  When  xcV  (47)  does  not  simply 
in  general.  It  describes  an  anisotropic,  inhomogeneous,  spa¬ 
tially  dispersion  medium. 

Let  us  examine  how  (53)  reduces  to  some  more  familiar 
expressions  in  some  special  situations.  We  will  assume  that 
V  is  infinite  through  the  remainder  of  this  section.  First 
we  will  consider  the  case  when  the  density  is  constant,  i.e., 
p(s_)=p.  In  this  case  the  permittivity  is  translationally  in¬ 
variant  or  homogeneous.  To  see  this  we  substitute  (21)  into 
( 53)  and  we  perform  the  integrations  over  s  and  k ' .  We  ob¬ 
tain 
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(54) 


c  =  I  + 

kt 


^dx'  J' dk 

when  ^(k,k)=  f  P_  (<*>_)  t  (k  ,k  '  ;w_) 

**  A  rr  _ 


eik-  <x-x’>  £(krk, 


(55) 


Since  the  integrand  is  a  function  of  x-x'  the  permittivity-  is 
translationally  invariant  however  it  is  still  anisotropic  and 
spatially  dispersive. 

Another  special  case  of  interest  is  when  t  is  scalar, 
i.e.  £(x,x')  =  t ( x , x ' ) I .  This  occurs  when  the  scatterers  are 
spherical.  Then  the  permittivity  is  isotropic  but  inhomo¬ 
geneous  and  spatially  dispersive. 

The  last  special  case  to  be  treated  is  when  the  wave¬ 
length  is  large  compared  to  the  size  of  a  scatterer.  Here 
the  particle  can  be  treated  as  an  electric  dipole.  Its  equi¬ 
valent  source  distribution  is  given  by 

-Seq  =  ^o^eq  =  iaiyo  ( -i^£<S  ( x)  )  (56) 

where  p  is  the  electric  dipole  ncment  of  the  scatterer  and  6  (x)  is  tne 
Dirac  delta  function.  The  dipole  moment  is  related  to  the  incident  field  eQ  by 
the  polarizability  tensor  a  [Jones,  196^] 

P  “  e0  g*£o  (57) 

Using  (57)  in  (56)  and  comparing  it  with  (14) ,  we  find 

T  =  k^gfi (x)  (58) 

or 

t (x, x ' )  =  kQa6(x)6(x')  (59) 

Now  putting  (59)  in  (19) ,  we  obtain 

£  (k  ,  k  '  )  -  k^a/ (  2tt  )  3  (60) 
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Thus  we  see  in  the  dipole  limit  t  is  independent  of  k  and  k_‘  . 

Since  we  have  an  expression  for  t  in  the  low  frequency 
or  dipole  limit,  the  special  form  of  e  can  be  easily  obtained 
Using  (5  9)  in  (54)  ,  we  find 

|  =  I  +  p  (x)  a  (61) 

Thus  in  the  low  frequency  case,  the  permittivity  is  no  longer 
spatially  dispersive  however  it  is  still  anisotropic  and  in¬ 
homogeneous  , 


CORRELATION 


In  this  section  we  will  calculate  the  correlation  of  the 

electric  field.  Rather  than  following  procedures  used  to  find 

the  coherent  wave,  the  distorted  Born  approximation  will  be 

employed.  This  is  a  single  scattering  approximation  where 

the  scatterers  are  assumed  to  be  embedded  in  the  equivalent 

medium  which  has  been  found  in  the  previous  section.  The 

method  is  useful  when  the  fractional  volume  is  small  (NV  /V<<1) 

P 

and  the  albedo  of  a  single  particle  is  small.  The  later  con¬ 
dition  implies  that  the  energy  absorbed  by  a  particle  must  be 
much  larger  than  the  energy  scattered  by  it. 

We  start  by  considering  a  volume  V  of  equivalent  medium 
surrounded  by  free  space.  There  are  N  particles  embedded  in 
V  as  shown  in  Figure  1.  The  scattered  field  due  to  the  ith 
particle  can  be  calculated  by  modifing  (39).  We  assume  that 
the  incident  field  on  the  particle  is  the  mean  field  <E>  and 
that  the  free  space  operator  i  is  replaced  by  the  equivalent 
medium  operator  as  given  in  (48) .  We  have 

N  .  .  .  N  , 

E  *  E  E  1  ■  Z  fc  -T.*<E>  .  (62) 

S  i=l  S  i=l  "  ==1 

Before  proceeding,  we  point  out  that  our  main  interest 
in  finding  the  correlation  of  the  field  is  to  use  it  to  cal¬ 
culate  the  backscatter ing  cross  section.  Since  this  cross 
section  is  related  to  the  correlation  of  the  field  fluctua¬ 
tions,  we  now  define 

Ef  «  E  -  <Eg>  ,  <Ef>  *  0  (6  3) 

Now  computing  the  correlation  of  the  fluctuating  field,  we 


obtain 


mum  a  m 


E,(x)E*(x)>  =  ^  E  ( x )  E  *  ( x )  >  -  <  E  (x)><E*(x)> 
— r  —  — f  —  — s  —  — s  —  — s  —  s  — 


'64  ) 


where  z*  is  the  conju_.  te  of  z.  Putting  (62)  in  (64)  and  no¬ 
ting  that  a  portion  of  <E  £*>  cancels  with  <E  ><E*>  if  we  use 

— s  —s  — s  — s 

the  fact  that  N>>1.  We  find 


<Ef  (x)  E*  (x)  >  =  f  <%  (&)  <Ef  (x)  E£  (x)  >[j 

-^4  7r  — 


(65) 


where 


A  /*  /v 

Ef(x)E^(x)>^  -  J  ds  p(s)  £(x,s)£*(x(s)  (66^ 


with 


>-l 


(67) 


£(x,s)  =  t  ’T(s)  «< E> 
here  we  have  separated  the  average  into  rotation  and  coordi¬ 
nate  space  averages,  thus  introducing  the  conditional  expecta¬ 
tion,  <E^E*>  ,  with  respect  to  given  to. 

~£—:t  m  ~ 

To  write  (67)  more  explicitly,  we  introduce  the  dyadic 
Green's  function  g(x,x')  for  the  operator  It  satisfies 

£-£(x,x')  =  Io(x-x')  (68) 

+  G  -  outgoing  as  ix.!-*00 
where  is  given  in  (4  8)  .  Now  (6  7)  becomes 

S(x,s)  =  _/dx'g(x,x’)  •  ./dx"|(x,-s,V’"S)  -<E  (x")  >  (69) 

The  expression  simplifies  greatly  in  the  low  frequency  limit. 
Assuming  that  t  is  given  by  (59)  and  using  this  in  (69)  gives 


Six,  s)  =  -^0S(x,s)  *g-<E  (s)  > 


(70) 


BACKSCATTERING  COEFFICIENTS  FOR  A  HALF  SPACE  OF  DIPOLES 


To  illustrate  the  application  of  the  methods  developed 
in  the  previous  sections,  we  will  calculate  the  backscattering 
coefficients  from  a  half  space  of  scatterers  that  are  small 

compared  to  wavelength.  We  will  alsp  assume  that  the  density 

/ 

of  scatterers  p  is  constant.  The  physical  configuration  is 
shown  in  Figure  2 .  There ,  we  have  shown  the  direction  of  the 
incident  wave  and  the  polarization  vectors  h°  and  v°  repre¬ 
senting  horizontal  and  vertical  polarizations  respectively. 

To  compute  the  scattered  field  using  the  distorted  Born 
approximation,  we  must  first  calculate  the  mean  field  in  the 
half  space  containing  the  particles.  In  the  low  frequency 
approximation  the  mean  wave  is  computed  by  replacing  the  par¬ 
ticles  with  an  equivalent  medium  having  relative  permittivity 
tensor  £~I+p  g  and  free  space  permeability  .  The  usual  con¬ 
tinuity  conditions  associated  with  macroscopic  Maxwell's 
equations  are  assumed  to  hold  at  the  interface  z~0 . 

Before  proceeding  we  would  like  to  emphasize  that  the 
scatterers  are  sparsely  distributed  or  that  the  fractional 
volume  they  occupy  is  small  (NV  /V<<1)  -  a  condition  necessa- 
sary  for  the  validity  of  the  mean  equation.  This  restriction 


is  reflected  in  the  equivalent  permittivity  tensor.  Small 
fractional  volume  requires  that  (pa.  .  |<<1  where  a .  .  are  tne 
components  of  3.  We  can  exhibit  the  dependence  of  this  con¬ 


dition  on  the  fractional  volume  6=NVp/V~pV  explicitly  hy  in¬ 
troducing  a  normalized  polarizability  tensor  a  as  follows: 


i  -  s/vp 


•  •  .»,», .. ......  ^jJtwhik.nVi  ....Am-  *'*|iii 


Figure  2  Incidence  wave  on  half  space  of  unif  ox-mily 
distributed  particles 
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where  we  can  show  that  the  components  of  a  remain  bounded  as 
V  -*■().  Now  the  permittivity  can  be  written  as 

|  =  I  +  5g  (72) 

and  thus  we  have  a  small  parameter  for  ordering  purposes. 

Although  we  are  able  to  carry  out  the  calcualtion  of  the 
mean  wave  for  an  arbitrary  average  polarizability  tensor  § , 
it  is  convenient  to  partially  specify  the  angular  probability 
density,  ( oj )  in  order  to  make  g  diagonal.  First  we  choose 
a  spherical  coordinate  system  of  mutually  orthogonal  unit 
vectors  r ^ ,  0^  and  .  The  position  of  these  vectors  is  com¬ 
pletely  determined  by  the  spherical  angles  0  and  <f>  as  shown 
in  Figure  3.  Now  ws  align  the  principal  axes  of  the  scatter- 
er  along  these  unit  vectors.  Then  we  write 

3  =  a^.  r°r°  +  aQ0°8°  +  a  4>°<p°  (73) 

=  r  —  --  0—  —  <p—  — 

By  using  the  usual  transformation  between  spherical  and  car¬ 
tesian  coordinates,  (73)  becomes 

3  3 

3  =  E  E  a  ( 0  ,  <f>)  x?x?  (74) 

1=1  3—1  i  3  J 

where  x.=x,  x2=y  and  x3~z  and  'X-°  and  r.°  dre  cartesian  unit 

vectors.  The  relationship  between  the  a  x  and.  ar  ,ag,a,  are  giv6n  in  Appendi 

i  j 

Now  assuming  that  the  random  variables  0^  and  4^  are  in¬ 
dependent,  we  write 

pn(«)  *  P0(6)p*(4)  (75) 

Averaging  (74),  we  have 


3  3 

Z  Z 
i“l  j=l 


X?X° 

-x-3 


where 
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<E  (x)  >  =  A  ex— (86) 
where  <_  =  <_t  +  sc  ^z°  .  In  order  to  match  fields  at  the  inter¬ 
face,  the  transverse  phase  velocity  of  the  incident  and  trans¬ 
mitted  waves  must  be  the  same,  thus  <_t=kt=k0sin60x°  .  Putting 
(06)  in  (85),  we  have 


2  2— 

k_x  (<_xA)  +  k^A  +  6kc>a*A  =  0 


(87) 


Representing  A  in  cartesian  components  (87)  can  be  written  as 


where 


xx 


r  2  .2 

z  0  xx 

0 

kosinVz 

V 
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"kne 
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-1 

—  — 

1  +  6a  , 

XX 

6 

yy 

2a 

=  cos  0Q 

+  6a  , 

yy 

Bzz  = 

=  0  (88) 


0  22 
(8  9) 


Since  (88)  is  a  homogeneous  system,  the  determinant  of  coef¬ 
ficients  must  be  zero  for  a  solution  to  exist.  This  condition 

determines  the  allowable  values  of  k  .  We  find 

z 


< 


z 


+kn/e 

-  0  yy 


(90  ) 


K  =  +  <v>  =  +k 
Kz  -  z  -  0 


Vii_ 

sin^e0+Bzz 


1/2 


(91  ) 


where  the  superscripts  h  and  v  have  been  used  to  designate  the 
propagation  constants  associated  with  horizontal  and  vertical 
polarizations.  For  an  incident  wave  that  is  not  grazing,  i.e. 
6 q /tt / 2 ,  expressions  (90)  and  (91)  can  be  simplified  using  the 


small  6  parameter.  We  have 
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+  0(52) 


(92) 


(h) 


=  kQ (cose„  + 


6a 

vv 

0  2 cos  S  , 


<  (V)  =  kn 

z  0 


s(a2Z+COs2e0<aXx‘aZ2! 

cos60  +  i\ - ssse: - 


+  0  (c  )  (93) 


We  cannot  calculate  these  propagation  constants  to  higher  ac¬ 
curacies  than  0(6)  since  the  original  mean  equation  has  only 
been  found  to  this  accuracy.  We  note  that  since  we  are  con¬ 
sidering  lossy  particles!,  the  k  are  complex  and  thus  the 
mean  wave  will  decay  away  from  the  interface. 

Next  we  calculate  the  amplitude  coefficients  for  the 
mean  wave  of  both  polarizations.  We  have 


<  E (x,h) >  = 


(e 


ik-z 
8+r.  e 


-ik. 


ikt. 


x 


)  e 


•  (h ) 

K  ’X 


and 


<E  (x , v)  >  = 


V 


ik,. 

.  lkzz  .  -lkzz, 

Ke  0+rve  0  )e 


iK  -x 
(A  x°+A  z  0 ) e  —  — 

x—  z— 


z<  0 


z>0 


(94) 


z<  0 


z>0 


(95) 


where  we  have  introduced  a  reflected  wave  in  the  free  space 
medium  at  the  specular  angle.  Now  by  using  the  fact  that  the 
tangential  <E>  and  <H>  must  be  continuous  at  the  interface, 
the  unknown  reflection  and  transmission  coefficients  can  be 
calculated.  Since  the  major  effect  of  the  equivalent  medium 
is  to  produce  exponential  decay,  we  next  expand  the  coeffi¬ 
cient  for  small  6  and  we  keep  only  zero  order  terms.  We  find 

that  r  -0(a)  and  thus  it  can  be  neglected.  The  transmitted 
<3 


mean  fields  are: 

i  (k  2  z+kt  ’2£) 

<E(x,q}>  =  q°e  0  +  0(6)  ,  z>0  (96) 

qe  {h,v} 

Proceeding  with  our  developemnt,  we  now  relate  the  transverse 
Fourier  transform  of  the  correlation  to  the  backscattering 
coefficients.  This  is  done  within  the  context  of  the  distor¬ 
ted  Born  approximation  developed  in  the  previous  section.  We 
start  by  taking  the  transverse  Fourier  transform  with  respect 


to  x  and  x  of  (65) ,  (66) 

and  (7  0  )  .  We  have 

<E£(kt,Z)E|.(kt,S)>  =  ^ 

A  'x  A  A 

duj_  pr2(w)<Ef  (kt,z)E*(kt  ,z)  >w 

.TT  —  — 

(97) 

•*-  A  A 

<E^(k,z)Ef(k,z)  >  *  p 

—X  —  u  —1  — t  UJ_ 

/•  -v  **  vs  /v 

J  ds  £(kt,z,s) -£*(kt,z,s) 

(98) 

where 

£(k_t,z,s)  = 

k^G  (kt ,  z,s)  (s)  > 

(99) 

Since  we  will  only  require  Q  when  x  and  s  are  in  the  equivalent 
medium  and  since  the  reflection  at  the  interface  is  small  we 
can  replace  G  by  the  dyadic  Green's  function  for  an  infinite 
equivalent  medium,  i.e., 

£(x,s)  =  S(oo)  (x-s)  +  0(6)  (100) 

We  have  written  the  infinite  space  Green's  function  in  terms 

of  x-s  since  it  is  translationally  invariant.  We  then  have 

~  (ca)  ~  (*>)  -ik  ■  £ 

S  °°  =  5  *  (kt,z-s)e  (101) 

Now  by  putting  (96),  (99),  (100) and  (101) in  (98),  by  Integra- 

.a 

ting  over  3.t  and  by  setting  z=z~0,  we  have 

<Ef  (kt»0)E*(kt,0)  =  S(kt,q|_m)  5  (kt-kt)  (102) 


where 


S(kt,q|£)  =  (  2tt)  2ok^V^  J"  ds  (Q  ( (kfe  > -s )  '  a  ’  g°  )  (g  ^  (kt'-s) 


-2  lime  (q)  s 


*  a  *  •  £  °)  .  e 


(103) 


Here  S(kt,qju)  is  the  transverse  dyadic  spectral  density  at 
the  interface  assuming  u  is  fixed.  The  normalized  polariza¬ 
bility  has  been  introduced  by  using  (71). 

By  using  the  results  of  Appendix  C  and  by  noting  that 
§  ( k  , q | we  obtain  the  backscattering  coefficients 


kocos20o 


£°*g(-k  ,q\u)‘P°  ,  p#qe  {h,v] 

^0 


(104) 


To  evaluate  the  integral  of  5  in  (103) ,  we  will  need  the 
transformed  Green's  function.  To  obtain  it,  we  first  write 


the  governing  equation  for  S 


It  is 


[Vx(Vx|)  -  k^J+Sa)]  ‘  G(W)  (x)  =  |6(x)  , 


(105) 


Using 


' (x)  -  outgoing  as  |x 


5 (0o)  (x)  = 


(106  ) 


in  (10  5)  ,  we  find 


[(jcxjcxl)  +  kQ  ( I+6|)]*  g  (<_)  =  -I 


(107  ) 


Then 


^  (  OO) 

s  ^t'z! 


d<2  gtOe1-’- 


(108  ) 


To  simplify  the  remaining  computation  for  Q  we  note  from 

(  CD  ) 

(104)  and  (103)  that  Q  will  only  be  required  for  £=-k  = 

0 

-kQsine  q x° . 

Inverting  (107)  and  performing  the  integral  in  (107)  by 


the  method  of  residues,  we  obtain 


("-t0' 2)  = 


sin60tc^  a(z) 

0  +  - %r- - (x°z0+z0x°) 

0 


( Jv)  k2g  . 

Uz  kOSxx) 


ei^2V)|Zl  i<Zh)i2 

2i^  "  + 

z  z 


fv°  (109) 


where 


a  (z) 


-  {.; 


(110) 


and  the  B's  are  defined  in  (89).  Approximating  the  coeffi¬ 
cients  to  zeroth  order  in  6,  we  have  the  simplified  expression 


i<2V)  l2l  ilC2h)  lzi 
S<“>  (-kt(),2)  =  rr fur-cos^  + 


(111) 


If  we  use  this  in  (103) ,  perform  the  integration  and  use 
the  result  in  (10 4) ,  we  have  our  final  form  for  the  backscat- 
tering  coefficients.  It  is 


6k*V  | a  T 

r°  «  .°.  .B<l! _ 

^  8ir(ImK^^  +  Ira<^) 


(112) 


where 


*  /  d“i 


£e  *a •£0 


(113) 


and  K  ,  se  (h,v}  are  given  in  (92)  and  (93).  The  dependence 
z 

of  | a  |  on  angle  of  incidence  is  worked  out  explicitly  in 
Appendix  B  for  scatferers  that  are  distributed  uniformly  in 
the  <f>  coordinate. 


The  final  result  given  in  (112)  can  be  expressed  in  terms 


of  scattering  cross  sections  of  individual  particles.  Using 
(28) ,  (60)  and  (71) ,  we  have 


f  =  p0,f*q°  -  2tt^  S  =  k^V  a  /4tr 
pq  *  =  M  pq  0  p  pq' 


(114) 


Now  by  recalling  that  the  backscattering  cross  section  from 

a  particle  a^^=4ir!f  j^,  (112)  becomes 

*  pq  ‘  pq1 


pc<*> 

■o  =  _  pq 


Pq  2ImK:(P)+2Im^(q) 


p,qe  {h ,v} 


(115^ 


z  z 

Following  Attema  and  Ulaby  [1978]  we  can  give  a  one  dimen¬ 
sional  interpretation  of  (115)  .  If  we  rewrite  (115)  as 


pq 


-  f°  dz  pa.(b>e 


-2Im< 


(q) 


-2Imic 


(P) 


pq 


(116) 


we  can  view  the  scattering  as  being  decomposed  into  scattering 
from  slabs  of  width  dz.  An  intensity  of  exp(-2lm<^q^  |z|)  is 
incident  on  the  slab  located  at  z.  The  incident  intensity  is 

/  U  \ 

backscattered  with  reflectivity  factor  per  .  The  backscat- 

pq 

rered  wave  then  decays  as  exp (-2Imic^q^  |z|)  until  it  reaches 
the  interface. 
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DISCUSSION  AND  NUMERICAL  EVALUATION 


In  this  section,  we  will  first  discuss  several  general 

properties  of  that  are  independent  of  the  particular  scat- 

terer  chosen.  Following  this  discussion,  we  use  our  method 

to  model  a  forest  canopy  by  a  collection  of  lossy  dielectric 

discs.  The  theoretical  curves  computed  from  this  model  are 

then  compared  with  some  experimental  data. 

Because  of  the  simple  dependence  of  a“g  on  the  medium 

properties  and  incidence  angle  certain  general  observations 

can  be  made  that  are  independent  of  the  particular  nature  of 

the  scatterer.  First,  we  note  that  c°^  as  given  by  (115)  is 

independent  of  the  density  of  scatterers  p.  This  follows  di- 

( s ) 

rectly  from  (92)  and  (93)  where  we  see  that  the  im«r  ,  se  (h,v) 

are  directly  proportional  to  p.  Thus  the  linear  p  dependence 

in  the  numerator  of  (115)  is  cancelled  out  by  the  denominator. 

Second,  we  note  that  o°  *  cos8..  This  is  the  same  result  as 

nn  u 

predicted  by  the  scalar  theory.  Finally,  we  note  that  o'^=°^r 
at  normal  incidence  (0°=*O).  This  is  an  expected  result.  Since 
the  scatterers  are  uniformly  distriouted  in  ,  the  two  polari¬ 
zations  see  the  same  medium  at  normal  incidence. 

We  now  proceed  to  model  a  forest  canopy  by  a  collection 
of  leaves.  The  leaves  are  in  turn  assumed  to  be  lossy  dielec¬ 
tric  discs  as  mentioned  previously.  The  discs  have  radius  a 
and  thickness  h.  Typical  dimensions  are  radii  of  one  to  sev¬ 
eral  centimeters  and  thicknesses  of  tenths  of  a  millimeter. 

The  electrical  properties  of  discs  can  be  characterized  by 
their  normalized  polarizability  tensor  a  when  the  wavelength 
is  large  compared  to  the  disc.  From  Jones  [1964]  and 


Van  de  Hulst  [1957]  the  polarizability  of  a  disc  along  its 
principal  axes  is  given  by 


V  1+A 
P 


(117) 


when  r,6f<p  are  defined  in  Figure  3. 

Because  of  the  large  volume  of  water  present  in  vegeta¬ 
tion,  we  can  usually  assume  |e  |>>1  in  the  microwave  region. 
Using  this  assumption  in  (117),  we  find  that  j aQ | = | a^ | >> | j . 
This  inequality  can  now  be  used  to  simplify  the  scattering 
cross  section  of  (115).  We  find  that  a®  =c|  j  2v^/Imer  • 

Thus  it  follows  that  the  magnitude  of  the  backscattering  cross 
sections  are  directly  related  to  the  volume  and  complex  die¬ 
lectric  constant  of  the  discs  in  a  simple  manner.  Therefore, 
as  leaves  grow  and  as  their  moisture  content  changes  these 

effects  should  be  observable  by  measuring  o°  at  different 

pq 

periods  of  the  growing  season. 

Before  computing  the  backscattering  cross  sections  as  a 
function  of  incidence  angle,  we  will  require  the  relative  di¬ 
electric  constant  of  the  leaves  and  the  angular  distribution 
of  leaves.  First  the  relative  dielectric  constant  is  considered 
Our  calculation  follows  that  of  Fung  and  Ulaby  [1978]  who  in  turn  have 
based  his  results  upon  de  Loor  [1968]  and  Carlson  [1967], 

They  model  the  leaves  as  a  mixture  of  water  and  solid  materi¬ 
als.  For  illustrative  purposes  we  have  chosen  50%  water  and 
50%  solid  for  our  calculations.  By  using  (3)  and  (4)  of 
Fung  and  Ulaby  [1978]  at  a  frequency  of  1.1  GHz  we  find  that  er= 
30.8+H.8.  our  choice  of  frequency  has  been  motiv  ated  by  experimen¬ 
tal  results  that  appear  in  the  literature.  We  have  chosen  to 


compare  our  results  with  those  of  Bush,  et  al  [1976]  who  has 


measured  a®,  p,qe  {h,v}  from  forests  for  frequencies  1-18  GHz. 

ir  Si 

Because  of  the  dipole  approximation  made  in  our  model,  only 
the  lowest  frequency  (l.lGHz)  Ulafcy  measured  was  used  for 
comparison  purposes. 

The  angular  distribution  of  leaves  will  now  be  considered. 
Field  measurements  of  leave  orientations  have  been  made  by 
Smith  [1973]  and  others.  It  has  been  found  that  the  leaves 
are  distributed  uniformly  in  the  4>  coordinate  (Figure  2) .  The 
distribution  of  leaves  with  respect  to  0  is  more  vegetation 
type  dependent.  Several  are  given  by  Smith  [1977].  Since 
no  measurements  of  this  type  exist  for  the  Ulaby  data,  we 
have  assumed  that  $  is  uniformly  distributed.  For  0  we  have 
considered  the  following  two  density  functions: 


°<  6<A6^ 
A  0  /f  <_u<t 


(118) 


or 


P0(0)  -  ( 


j  2A0X  ' 


TI  -A0.J  <  0<  TT  +  A0 


elsewhere 


(11?) 


In  (11 8)  when  A 8  is  small,  the  lec.ves  are  approximately 
parallel  to  the  interface  (z=0) ;  when  ,  they  are  un¬ 

iformly  distributed  in  0.  In  (119)  when  A0j,  is  small,  the 
leaves  are  perpendicular  to  the  interface;  when  A0^  is'  tt/2 
they  are  uniformly  distributed  in  0. 

The  numerical  calcualtionc  are  presented  in  Figures  4-9. 
In  these  figures  the  backscattering  coefficient  is  plotted  as 
a  function  of  the  angle  of  incidence  9_.  In  Figures  4-7  we 


ANGLE  OF  INCIDENCE  (DEGREES) 


Figure  10 


have  used  the  angular  distribution  given  in  (118)  while  in 

Figures  7-9  we  have  used  (119).  Figure  7  corresponds  to  a 

uniform  distribution  over  all  9  and  thus,'  for  this  case, 

(118)  and  (119)  give  the  same  results.  One  should  note  that 

since  |e  |>>1,  a  change  in  f,  V  or  £  just  shifts  the  level 
r  p  r 

of  the  curves  but  it  does  not  change  their  shape.  Their  shape 
only  depends  on  Pg(0). 

The  following  trends  are  observed  in  Figures  4-7:  First, 


is  always  greater  than  a' 


Their  difference  increases 


as  A  6^  becomes  smaller.  Second,  the  cross  polarized  back- 

scatter  becomes  smaller  as  A becomes  smaller.  In  Figure 

7-9,  we  observe  that:  First,  becomes  greater  than  ahh  33 

A 9^  is  decreased;  second,  the  curve  for  o£v  tends  to  flatten 

out  as  A9^  is  increased;  and  third,  the  difference  CTj1j1“°hv  at 

9q=0  becomes  smaller  as  A9j  increases. 

A  comparison  of  our  theory  with  the  experimental  results 

of  Bush,  et.al  [1976]  is  made  in  Figure  10.  There,  we  have 

plotted  our  Figure  8  along  with  his  data  for  Kansas  deciduous 

trees  measured  in  the  springtime  at  a  frequency  of  1.1  GHz. 

Figure  8  was  chosen  since  it  most  clearly  appears  to  follow 

the  trends  of  the  data  i.e.,  flat  cross  polarization  and  a0  > 

rrv 


Although  our  theory  follows  the  trends  of  the  data,  it 
is  clear  from  the  results  that  additions  to  our  model  .should 
be  made.  The  fact  that  a-^XJhh  mostL'  likely  due  to  the 
vertically  oriented  tree  branches  other  than  due  to  the  leaves 
that  tend  to  be  parallel  to  the  interface.  In  addition,  an 
examination  of  the  numerical  results  shows  that  the  skin  depth 


for  the  mean  wave  is  large.  Thus  at  a  frequency  of  1  GHz  the 
under lyir.9  ground  should  be  taken  into  account. 


CONCLUSIONS  AND  RECOMMENDATIONS 


We  have  developed  and  analyzed  a  discrete  scattering 
model  for  vegetation.  The  model  replaces  the  vegetation  by 
simple  scattering  objects  such  as  discs,  spheres  or  rods  so 
that  leaves,  branches  and  trunks  can  be  modeled.  The  objects 
are  given  random  placement  and  orientation  in  order  that  the 
complicated  features  of  the  individual  scattering  objects  are 
averaged  out  and  simple  expressions  are  developed  for  the 
backscattering  coefficients.  The  model  directly  relates  the 
backscattering  coefficients  to  easily  measured  medium  param¬ 
eters  such  as  leaf  size,  density  and  moisture  content.  The 
theory  as  presented  also  takes  into  account  depolarization 
effects  that  agree  with  experimental  data. 

At  present  only  the  leaves  in  the  forest  have  been  in¬ 
cluded  in  this  model  and  no  effects  of  the  ground  have  been 
taken  into  account.  We  recommend:  first,  that  the  effect 
of  the  ground  be  included.  If  this  is  done  agricultural  crops 
can  be  modeled.  Second,  the  numerical  calculations  indicated 
the  branches  in  forests  will  make  an  important  contribution 
of  microwave  frequencies.  They  can  be  modeled  by  dielectric 
cylinders  and  averaged  over  position. 

Finally  we  feel  that  the  inverse  problems  related  to  the 
above  approach  should  be  investigated.  Questions  such  as: 

Can  the  probability  density  function  of  leaf  inclination  angles 
be  determined  from  backscatter  data?  Once  the  distribution 
of  inclination  angles  is  known,  forest  or  crop  identification 
will  be  much  easier. 
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APPENDIX  A 


Relationship  Between  the  Transition  Operator  and  the  Scatter¬ 
ing  Amplitude 

To  find  the  relationship  between  f  and  t,  we  start  with 
(15) .  Using  (23)  and  (26)  ,  we  can  write  (15)  in  terms  of 
dyadic  incident  and  scattered  wave 

(1A) 


e  =  L  1 *T • en 
=s  -  =  =0 


Next  we  use  (16)  in  (1A)  along  with  the  free  space  dyadic 
Green’s  function  £  for  L  \  We  have 

e  (x,_i)  =  Jdx'L  (x,x' )  'fdx"  t  (x'  ,x")  *§Q  (x"  ,i)  (2A) 

Here  the  free  space  dyadic  Green's  function  is  given  by 

ik0 1*“*' I 

-  (1+  fr’wiwr  <3A) 

*0 

To  obtain  _e  in  the  radiation  zone,  the  far  field  expression 
for  £  will  be  required.  It  is  [Twersky,  1967] 
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Now  putting  (4A)  and  (25)  in  (2A) 
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Finally  employing  (19)  in  (5A)  ,  and  comparing  w:tth  (27)  we 
obtain  the  required  result 


APPENDIX  B  ~  Polarizability  Statistics 


In  this  appendix  calculate  the  mean  square  polarization 
statistics  used  in  (112)  .  The  calcualtion  will  be  performed, 
with  the  assumption  that  the  scatterers  are  uniformly  distri¬ 
buted  in  $. 

First  we  derive  the  components  of  the  polarizability 
tensor  in  cartesian  coordinates  in  terms  of  the  principal  axis 
components.  The  unit  vectors  r1 0 ,  _9°  and  j).0  are  related  to  x°, 
and  z°  as  follows: 

r°  =  sinQcosdx0  +  sin0sin<f>y°  +  cosSz0 

0°  =  cos9cos$x°  +  cos0sin<)>y.0  -  sin0£°  (IB) 

jp°  =  -sin$x°  +  cos<j>y° 

Using  this  in  (73),  the  cartesian  components  of  (74)  are: 

2  2  2  2 
a  „  =  (a  sin  0  +  a0cos  0)cos  <p  +  a,  sin  $ 
xx  i  y  <p 

2  2  2  2 

a  =  (a  sin  0  +  aQcos  0)sin  $  +a,cos  <p 
yy  r  y  <p 

2  2 

a  ■=  a  cos  -0  +  aQsin  9 

zz  r  9 

(2b; 

2  2 

a  „  =  [(a  sin  9  +  aQcos  9)  -  a ,  ]  cos<}>sind 
xy  r  y  <p 

a  *  (a  -a fl)  sin0cos9coscf> 

X  2  1  u 

a  -  (a  -aQ)  sin9cos9sin4) 
yz  r  a 

The  other  components  are  gotton  from  the  fact  that  a  is  a 
symmetric  dyadic. 

New  we  will  obtain  the  components  of  a  in  the  polariza¬ 
tion  directions  in  terms  of  the  cartesian  components.  We 
have  using  (8  4) 
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J-.-7  ■  r-rr  *  cos2«oaXx  -  2oos9osineoaxZ  +  sin2Vz2 

Following  this,  the  mean  square  polarizabilities  are  com¬ 
puted.  When  we  assume  ag=a^=:afc  and  we  use  the  approximation 
that  a^>>a  ,  the  mean  square  polarizabilities  become 
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APPENDIX  C  -  Relationship  Between  the  Backscatter.ing 

Coefficients  and  the  Transverse  Spectral  Density 


We  start  by  considering  the  fluctuating  portion  of  the 
scattered  field,  ,  as  defined  in  (5?!  in  the  region  z  _>0 . 
This  field  can  be  viewed  as  arising  from  sources  on  the  in¬ 
terface.  So  chat  far  field  quantities  can  oe  found,  we  ini¬ 
tially  consider  that  portion  of  that  arises  from  sources 
contained  within  a  finite  region  A  on  the  interface.  The  ra¬ 
diated  field  from  the  region  A  will  be  denoted  by  E,  (x, q)  . 

A 

The  field  E^  (x,q)  can  be  related  to  the  interface  fluc- 
‘“A 

tuations  by  employing  a  plane  wave  expansion  in  the  region 
z<_0 .  We  have 

l  f  r  -ik  z+ik  (xt-x‘) 

E  -  ( x ,  q }  - 7/  dk  (  dx  lf(0,x/,q)e  ,  z<  0 

-*  ,lci 

where  Ef(0,xt,q)  is  the  fluctuating  field  on  the  interface 
due  to  an  incident  wave  of  polarization  q.  The  k  integral 
in  (1C)  can  be  asymptotically  evaluated  for  large  [x[  [Collin 
and  Sucker,  1969]  .  We  find 
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where  x  has  been  specialized  to  the  backseat ter  direction  9( 
and  k  is  given  by  (83) . 


The  backscattering  ccef f icients  are  now  defined.  They 


are 
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where  1.  (q)  is  the  incident  intensity  per  unit  area  with  po¬ 
larization  q,  I  (p,q)  is  the  average  intensity  at  the  obser¬ 
vation  point  with  polarization  p  due  to  an  incident  wave  with 
polarization  q.  In  view  of  (81),  0(q)=l,  We  also  have 

I  (p,q)  =<!  E  (X/q)*p.°!2>  (4C) 

S  rA 


We  can  complete  the  development  by  first  representing 
the  field  on  the  interface  by  its  Fourier  transform,  E^Oc^q): 

ik^  *  x. 


E  f  ( 0  ,  x. t ,  q )  = - 2 

(2t) 


y*Ef(k^,q)e  fcdkl 
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Now  by  using  (20  ,  (40  and  (50  in  (30  ,  we  have 
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Next  we  introduce  the  transverse  dyadic  spectral  density 
§  (k.t  »q)  whicn  is  given  by 


<^f  (llt'q)Ef  (h.t'q)  >  =  |(kt/q)  6  HO 

Note  that  this  definition  requires  that  the  process  Ef(kt,g) 
be  homogeneous  or  stationary  in  k ^ . 

By  using  (7C)  in  (60  and  carrying  out  the  integrations, 
we  obtain 
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where 


SD^kt»q)  =  £°  'I  (k_t  ,q)  -£c 
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To  cope  with  the  expandng  technology,  our  society  must 
he  assured  ol  a  continuing  supply  ot"  rigorously  trained 
and  educated  engineers.  The  School  ot  Engineering  and 
Applied  Science  is  completely  committed  to  this  ob¬ 
it  dive. 


